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Abstract 

The standard model (SM) plus a real gauge-singlet scalar field dubbed darkon (SM+D) is the simplest 
model possessing a weakly interacting massive particle (WIMP) dark-matter candidate. The upper limits 
for the WIMP-nucleon elastic cross-section as a function of WIMP mass from the recent XENON10 and 
CDMS II experiments rule out darkon mass ranges from 10 to (50, 70, 75) GeV for Higgs-boson masses of 
(120, 200, 350) GeV, respectively. This may exclude the possibility of the darkon providing an explanation 
for the gamma-ray excess observed in the EGRET data. We show that by extending the SM+D to a two- 
Higgs-doublet model plus a darkon the experimental constraints on the WIMP-nucleon interactions can be 
circumvented due to suppression occurring at some values of the product tana tan/3, with a being the 
neutral-Higgs mixing angle and tan/3 the ratio of vacuum expectation values of the Higgs doublets. We 
also comment on the implication of the darkon model for Higgs searches at the LHC. 
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I. INTRODUCTION 



The standard big-bang cosmology describes various astronomical and cosmological observations 
very well. To explain observational data, it requires two exotic unknown species beyond the standard 
model (SM) of particle physics, namely the dark matter and the dark energy, which in the energy 
budget of the Universe make up roughly 20% and 75% of the total energy density, respectively 
Although the evidence for dark matter has been established for many decades, the identity of its 
basic constituents has so far remained elusive. 

One of the popular candidates for dark matter (DM) is the weakly interacting massive particle 
(WIMP). Needless to say, the detection of a WIMP candidate is crucial not only for understanding 
the nature of DM, but also for distinguishing models of new physics beyond the SM. A variety of 
experiments have been and are being carried out to detect DM directly by looking for the recoil 
energy of nuclei caused by the elastic scatteringof a WIMP off a nucleon. Although there is some 
evidence for WIMPs from indirect DM searches [l| , there is as yet no confirmed signal for them from 
direct searches. Stringent bounds on the WIMP-nucleon elastic cross-section have been obtained 
from the null results of direct DM searches. The strictest limits to date on the cross section as 
a function of WIMP mass have been set by the recent XENON10 0] and CDMS II j^] experiments. 
They are shown in Fig. (TJ along with the expected sensitivities of a few future experiments. These 
upper bounds provide important restrictions on WIMP models. 

Among the many possible WIMPs, the lightest supersymmetric particle is the most studied 
candidate. Although this possibility has many attractive features, no direct experimental evidence 
for it has yet been discovered. It is therefore important to study and search for other possible 
WIMPs which can explain the DM relic density. The simplest model which has a WIMP candidate 
is the SM+D, which extends the SM by the addition of a real gauge-singlet scalar field D. This 
singlet field, which we will call darkon, can play the role of dark matter. The darkon as DM was 
first considered by Silveira and Zee j3| and further explored later by other groups 
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CDMS: 2004+2005 (reanalysis) +2008 Ge 
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FIG. 1: Current and projected upper limits for the spin-independent WIMP-nucleon elastic scattering 
cross-section as functions of WIMP mass. Plot courtesy of Ref. 
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In this paper we study the constraints on the SM+D and its two-Higgs-doublet extension 
(THDM+D) from the recent XENON10 and CDMSII experiments which set upper limits for 
the WIMP-nucleon spin-independent elastic scattering cross-section 0, S|- In the SM+D, we find 
that a darkon having a mass in the ranges of 10 GeV to (50, 70, 75) GeV for Higgs masses of 
(120, 200, 350) GeV, respectivejvis ruled out by the experimental bounds. This confirms some of 
the results obtained in Refs. |S|, |9j. Now, a WIMP candidate with a mass in the range from 50 to 
70 GeV can well explain the gamma-ray excess observed in the EGRET data [l(|. The darkon mass 
range ruled out by XENON10 and CDMS II would imply that the darkon model is not likely to offer 
an explanation for the EGRET excess. It is then interesting to see if there are possible avenues to 
evade the experimental restrictions on the WIMP-nucleon interactions and keep this darkon mass 
range viable. To this end, we consider two-Higgs-doublet models with a darkon and find that in 
the type-II THDM plus a darkon (THDMII+D) it is indeed possible for a darkon with a mass 
within the range of interest to satisfy the experimental requirements. This can happen because of 
the suppression of the darkon-nucleon interaction at certain values of the parameters of the model. 
We also find that a darkon having mass within this range can lead to considerable enhancement 
of Higgs decay widths through substantial invisible decays into darkon pairs. This could cause the 
Higgs branching fractions to SM particles to diminish and consequently affect Higgs signatures at 
colliders significantly. 

Before starting our analysis of the models, we would like to summarize the relic-density con- 
straints that any WIMP candidate has to satisfy. For a given interaction of the WIMP with SM 
particles, its annihilation rate into the latter and its relic density Q D can be calculated and are 
related to each other by the thermal dynamics of the Universe within the standard big-bang cos- 



mology 12[ . To a good approximation, 



H 107 ; x 109 */ - , x . H , (i) 

/<U™>pi (a ann i; rcl ) GeV ' ^xj 



where h is the Hubble constant in units of 100km/(s-Mpc), m pl = 1.22 x 10 19 GeV is the Planck 
mass, m D is the WIMP mass, Xj = m D /Tj with Tj being the freezing temperature, is the 
number of relativistic degrees of freedom with masses less than Tf, and (c ann v rel ) is the thermal 
average of the product of the annihilation cross-section of a pair of WIMPs into SM particles and 
the relative speed of the WIMP pair in their center-of-mass frame. 

The current Particle Data Group value for the DM density is Q D h 2 = 0.105 ± 0.008 13] . Using 



this number and Eq. ([TJ), one can restrict the ranges of Xf and (er ann f rc i) as functions of WIMP 
mass m D without knowing the explicit form of the SM-WIMP interaction. In Fig. [2] we show the 
values allowed by the 90% C.L, range 0.092 < Q D h 2 < 0.118. 



II. STANDARD MODEL WITH DARKON 



A. Brief description of the model 



Since the darkon field D in the SM+D must interact weakly with the SM matter fields to play the 
role of dark matter, the simplest way to introduce the darkon is to make it a real gauge singlet with 
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FIG. 2: Ranges of Xj and (a 
0.092 < n D h 2 < 0.118. 



as functions of WIMP mass corresponding to the 90% C.L. range 



a discrete Z2 symmetry so that it does not have SM gauge interactions and can only be created or 
annihilated in pairs. Requiring that the darkon interactions be renormalizable implies that D can 
only couple to the Higgs doublet field H . Beside the kinetic energy term ^d^D d^D, the general 
form of the other terms in the darkon Lagrangian is 0, [5J 

C D = -^D*-^D 2 -XD 2 H^H , (2) 

where A^,, m , and A are free parameters. Clearly, C D is invariant under the Z 2 transformation 
D — > —D, the SM fields being unchanged. The parameters in the potential should be chosen such 
that D does not develop a vacuum expectation value (vev) and the Z2 symmetry is not broken, 
which will ensure that the darkon does not mix with the Higgs field, avoiding possible fast decays 
into other SM particles. 

The Lagrangian in Eq. ([2]) can be rewritten to describe the interaction of the physical Higgs 
boson h with the darkon as 

C D = ~^-D D --D h - XvD h, (3) 

where v = 246 GeV is the vev of if, the second term contains the darkon mass m D = [rn^ + Xv 2 ) 
and the last term, —XvD 2 h, plays an important role in determining the relic density of the DM. At 
leading order, the relic density of the darkon results from the annihilation of a darkon pair into SM 
particles through Higgs exchange 0,0], namely DD — > h* — > X, where X indicates SM particles. 

Since the darkon is cold DM, its speed is nonrelativistic, and so a darkon pair has an invariant 
mass \fs ~ 2m D . With the SM+D Lagrangian determined, the /i-mediated annihilation cross- 
section of a darkon pair into SM particles is then given by j^j 

(j „ 1 = ^ £.r(^.) 

— ' (iml-mlf + rlml 2m D U 

where v rel = 2\p c ^\/m D is the relative speed of the DD pair in their center-of-mass (cm) frame, 
h is a virtual Higgs boson having the same couplings to other states as the physical h of mass m h , 
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but with an invariant mass y/s = 2m D , and h — >• Xi is any possible decay mode of h. For a given 
model, TliT[h — > XA is obtained by calculating the h width and then setting m h equal to 2m D . 

The darkon-Higgs coupling A for a given value of m D can now be inferred from the range of 
(& aim v rcl ) values allowed by the Q D h 2 constraint, as in Fig. PSJ once m h is specified. We show 
in Fig. [3] the allowed ranges of A corresponding to 10 GeV < m D < 100 GeV for representative 
values of the Higgs-boson mass. We note that it is possible for A to become larger than 1 when 
m D decreases, which would upset the applicability of perturbative calculation. Consequently, we 
display only the A < 1 regions in this figure. 




m D (GeV) 

FIG. 3: Darkon-Higgs coupling A in the SM+D as a function of darkon mass m D for Higgs mass values 
m h = 120, 200, 350 GeV. 



B. Effective Higgs-nucleon coupling 

The detection of dark matter on the Earth is through the recoil of nuclei when a darkon hits 
a nucleon target. This interaction occurs via the exchange of a Higgs boson between the darkon 
and the nucleon N in the t-channel process DN — > DN. Its cross-section then depends on not 
only A, but also the Higgs-nucleon coupling. Since the energy transferred in this elastic scattering 
is very small, of order 100 keV, one can employ a chiral-Lagrangian approach to obtain the effective 
Higgs-nucleon coupling. The Higgs-nucleon coupling has been studied previously 0, 0, and 
here we reorganize the results for convenience in evaluating possible cancelation among different 
contributions. 

The Higgs-nucleon coupling depends on the underlying Yukawa interactions of the Higgs boson 
with the quark degrees of freedom. The couplings of a Higgs boson T-L to quarks can be generically 
written as 

k 

L mH = -^^ m q ^' H - ( 5 ) 

i 

where the sum runs over the six quark flavors, q = u,d,s,c,b,t. In the SM k q = 1 for all g's, 
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whereas beyond the SM the A; 's may have other values. Specifically, k u = k c = k t and k d = k s = k b , 
the two sets being generally different, in the THDMII+D which we will consider later. 

Now, the effective coupling of % to a nucleon N = p or n has the form 

C-nnh = -QnnhNNH, (6) 

where Qnnh ls ^ ne Higgs-nucleon coupling constant. In the SM+D and THDMII+D, one then 
needs to evaluate the matrix element 

k k 
9nnh NN = {N\ — (m u uu + m c cc + m t tt) + —(m d dd + m s ss + m b bb)\N}. (7) 

We outline the derivation in the Appendix, and the result is 

where a nN is the so-called pion-nucleon sigma term, m N the nucleon mass, and m B the baryon mass 
in the chiral limit. As explained in the Appendix, m B is related to a^ N by Eq. flA3j) and numerically 
we adopt a Tjv = 45 MeV. It follows that 

g NNn ~ (l.217A; d + 0.493A; u )xl0- 3 . (9) 

Hence g^N^t can ^ e vanishingly small if k u ~ —2.47 k d . 



C. Darkon-nucleon elastic cross-section in SM+D 

The amplitude for the elastic scattering DN — > DN mediated by % is given by 

M cl ~ 2X9n ™ V NN (10) 
for the squared transfer momentum t -C m^. This leads to the cross section 

7r (m^, + m N ) 

the approximation (pd + PnY — ( m D + m 7v) 2 having been used. In the SM+D, the Higgs-nucleon 
coupling g^ffh is obtained from Eq. ([§]) by setting k u = k d = 1. Thus 

~ 1.71 x lO" 3 , (12) 

which is comparable to the values found in the literature [5], [l5| . 

Since A as a function of m D is constrained by the relic density, so is <r el . In Fig. H] we plot 
(J cl versus m D in the SM+D for A < 1 and representative values of the Higgs-boson mass. We 
remark that the dip of the m h = 120 GeV curve at m D = 60 GeV in this figure corresponds to the 
dip in Fig. [3] and the minimum of the denominator in Eq. (j3J) at m h = 2m D . When comparing 
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predicted elastic cross-sections with experimental limits, we assume in this paper that the local 
density of dark matter is fully accounted for by the darkon, and so no scaling is needed for the cross 
sections. In Fig. HJ we also plot 90% CL. upper limits on the WIMP-nucleon spin-independent 
elastic cross-section set by the XENON10 and CDMSII experiments 0, El]. One can easily see 
that the darkon mass ranges 10 GeV < m D < (50, 70, 75) GeV are ruled out for Higgs masses of 
(120, 200, 350) GeV, respectively. 

The analysis above suggests that the experimental constraints can be evaded if the Higgs-nucleon 
coupling gets sufficiently small due to (at least) partial cancelation among various contributions to 
it from the underlying Yukawa couplings. As we will see in the next section, this possibility can be 
realized in the two-Higgs-doublet extension of the SM+D. 
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FIG. 4: Darkon-nucleon elastic cross-section <t c1 in the SM+D as a function of darkon mass m D for Higgs 
mass values m h = 120,200,350 GeV, compared to 90%C.L. upper limits from XENON10 (black dashed- 
curve) and CDMS II (gray [brown] dashed-curve) . 




III. TWO-HIGGS-DOUBLET MODEL WITH DARKON 



A. Brief description of the model 



As the name implies, the THDM has two Higgs doublets H\^- Depending on how they couple 
to the fermions in the SM, there are different types of the THDM 16(. In the THDM I, only one of 
the Higgs doublets gives masses to all the fermions. In the THDM II, the up-type fermions get mass 
from only one of the Higgs doublets, say if 2, and the down- type fermions from the other doublet. 
In the THDM III, both H\ and H2 give masses to all the fermions. 

In the THDM I, since only one Higgs doublet generates the fermion masses, the Higgs couplings 
to fermions are those in the SM, up to an overall scaling factor. Therefore, the Higgs couplings in 
the THDM I+D are similar to those in the SM+D studied in the previous section and thus cannot 
help overcome the difficulties encountered. In the THDM III, there are flavor- changing Higgs-quark 
couplings which introduce too many parameters for the model to have predictive power. For these 
reasons, we will consider only the THDM II with the darkon field added (THDMII+D). 
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The Yukawa interactions of the Higgs fields in the THDM II are described by 16 

£ Y = -Q L \ U 2 H 2 U R - Q L \fH 1 V R - Lj&H^n + H.c. , (13) 

where Q, U, T>, L, and E represent the usual quark and lepton fields and A n,d,/ contain the Yukawa 
couplings. To effect the separate couplings of H 1 and H 2 to the down and up sectors, respectively, 
it is necessary to introduce a discrete Z' 2 symmetry, under which H 2 — > —H 2 and U R — > —U R , the 
other fields being unaltered. In terms of their components, the Higgs doublets are 

H„ = ±( ^ ), (14) 

where k = 1,2 and v k is the vev of H k . Here h k and I k are related to the physical Higgs bosons 
H + and A and the would-be Goldstone bosons w and z by 



h\\ ( cos — sin (3 \ ( w + \ ( I x \ / cos (3 —sm.f3\fz 

h\ J \ sin/3 cos/3 J \ H+ J ' V J 2 / ~ V sin/3 cos/3 J ^ 

with tan /3 = v 2 /v l , whereas h k can be expressed in terms of mass eigenstates H and h as 

/i 1 \ / cos a — sin a\ f H 



h 2 } \ sin a cos a I \ h 



(15) 



(16) 



Hence the angle a indicates the mixing of the two CP-even Higgs bosons. 

In analogy to Eq. (T5]) in the SM+D case, in the THDM II+D we have the renormalizable darkon 
Lagrangian 

C D = -^D*-^D*-(\ lH lH 1 + \ 2 H}H 2 )D*. (17) 

As in the SM+D, here we have again imposed the Z 2 symmetry under which D — > —D with all the 
other fields unchanged, and for the same reasons we need to keep it unbroken. This Lagrangian 
also respects the Z' 2 symmetry mentioned in the preceding paragraph. 

After electroweak symmetry breaking, Co contains the D mass and the DD(h,H) terms given 

by 

m 2 D = ml + (A-l cos 2 /3 + A 2 sin 2 f3)v 2 , 

C DDh = — (— A x sin a cos/3 + A 2 cos a sin. (3^)vD 2 h = —\ h vD 2 h, (18) 

C DDH = — (A x cos a cos/3 + A 2 sin a sm(3)v D 2 H = —X H vD 2 H, 

with v 2 = v\ + v 2 , but C D has no DDA term. Since m , A 1; and A 2 are all free parameters, we can 
treat the mass m D and the effective couplings X h H as new free parameters in this model. 

From Eq. ( TT3]) . the Yukawa interactions of h and H are described by 

- „„,.,, / cos a h sin a H \ - „ rr , ( sin a h cos a H 
C ffH = -U L M^U R + - V L M d V R 



sin (3 v sin /3 v J \ cos /3 v cos (3 v 

sin a h cos a H 

cos (3 v cos (3 v 



sin a h cos a H . 

- E L M E R - + - + H.c. . (19) 
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We have not written down the couplings of the CP-odd Higgs boson A to fermions because it does 
not couple to D and consequently A is irrelevant to our darkon relic-density study, not contributing 
to darkon annihilation. 

To evaluate the darkon annihilation rate, we also need the couplings of h and H to the weak 



bosons V = W ± , Z. From the kinetic sector of the THDM, we have 16 



i 2 777 ^ 777 ^ 

C VVH = * W + "W- + — Z»Z^ [h sin(/3 — a) + H cos(/3 - a)] . (20) 



v v 



B. Darkon-nucleon elastic cross-section in THDM II+D 



Using the formulas given above, one can express the cross-section of the darkon-nucleon elastic 
scattering in the THDM II+D as 



™2 2 / \ „THDM \ n THDM\ 2 

^ ~ ^ l \9 N Nh *h9nnh (21) 

-(m D + m N ) 2 \ m 2 h m 2 H J 



7T 



where, from Eq. ([8]), the nucleon coupling to Ti = h or H is 

THDM _ (uH 1,H\°tvN , uH m N , 4fc^ — 25fc^ m B 

9nnh ~ \ K u - K d ) + K d — + 27 ~ ' [ > 

The parameters are read off from Eq. ( fl9l) to be 

, , cos a , h sin a , „ sin a , „ cos a , , 

K u — ■ n ; K d — n i K u — ■ n 5 K d — o ■ K^ ) 

snip cosp snip cosp 

If both h and H contributed, the analysis would be complicated. For concreteness, in our 
numerical analysis we will neglect contributions from H by requiring X H to be very small or m H 
very large. Under this assumption, we can see that, since k\ d are free parameters, the situation 
in the THDM II+D can be very different from that in the SM+D. Here it is possible for the terms 
proportional to k!^ d in the Higgs- nucleon coupling gj^^ 1 to cancel. From Eq. f l2"2"j) , the cancelation 
condition is 

% = -tana tan/3 = 27^ + 8m B 

Since m B is related to cr nN by Eq. ( 1A3j) and since a nN is not well determined, with values within 



the range 35 MeV < a nN < 80 MeV having been quoted in the literature [15|, |l7|, |l8j , the value 
of k d /k^ has a sizable uncertainty. Nevertheless, we have checked that the main conclusion of this 
section below does not change for this a nN range. For defmiteness, in the following we employ 



a nN = 45 MeV [15J, as in Eq. Q9J), yielding 



f h ~ -0.405 . (25) 
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We have found that one can get the darkon-nucleon cross-section to be near the experimental 
bounds by allowing tan a tan /3 to deviate from the cancelation point, without violating the relic 
density constraint. For illustration, we choose tan a tan (3 = 0.45 and consider two sets of tana 
and tan /3 values satisfying this choice: (tana, tan (3) = (0.45, 1) and (tana, tan 0) = (0.45/30, 30). 
We note that, since low values of tan/3 are disfavored for low values of the charged-Higgs mass 
19[| , we will assume a large (> 250 GeV). We also note that tan /3 = 1 and 30 roughly 



span the range allowed by various experimental and theoretical constraints 19 



Since the relic density in this case is determined by the interaction of the darkon with h alone, 
the coupling X h can be extracted following the steps taken in Sec. [TTJ We show the resulting values 
of X h in Fig. [5j With determined, we can calculate the darkon-nucleon elastic cross-section, 
which is now given by 1 



m" N v 



\ „THDM 
A h ilNNh 



vr (m D + m N )' 



m 



(26) 



In Fig. |6]we plot a cl for our parameter choice above and A/ t < 1. We note that, as in the SM+D 
case, the dip of each of the m h = 120 GeV curves in this figure at m D = 60 GeV corresponds to the 
minimum of the denominator in Eq. (j4j) at m h = 2m D . In Fig. [6] we also plot the XENON10 and 



CDMS II upper limits, along with the expected sensitivities of a number of future experiments [11 



Evidently, there is parameter space in the THDMII+D that can escape the current experimental 
constraints. Future direct-search experiments can place more stringent constraints on the model. 



tanoc=0.45 and tan|3 = l 



tanoc=0.45/30 and tan|3=30 
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FIG. 5: Darkon-Higgs coupling X h in the THDMII+D as a function of darkon mass m D for Higgs mass 
values m h = 120, 200, 350 GeV and two cases with different tan a and tan /3 values satisfying tan a tan j3 = 
0.45. Only regions corresponding to X h < 1 have been plotted. 



It is worth remarking that, although the a and /? values chosen to make <?jv™ sufficiently small may at the same 
time cause gjfjpx to be enhanced, the contribution of H to a cl in Eq. (pH]) can still be neglected, as X H remains 
a free parameter of the model. As Eq. ([T5]) shows, X h H depend on the free parameters X l 2 , and so fixing does 
not imply fixing X H . Consequently, the product X H gj^f^ 1 can always be made as small as desired. In addition, 
as mentioned above, m H can be taken to be very large to suppress the second term in Eq. (|21[) further. 
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FIG. 6: Darkon-nucleon elastic cross-section a el in the THDMII+D as a function of darkon mass m D 
for Higgs mass values m h = 120, 200, 350 GeV and the two cases in Fig. [5l compared to the 90% C.L. 
upper limits from XENON10 (black dashed curve) and CDMSII (gray [brown] dashed curve), as well as 
projected sensitivities of SuperCDMS at Soudan (green dot-dashed curve), SuperCDMS at Snolab (gray 
[brown] dotted curve), and XENON100 (black dotted curve). 



IV. DISCUSSIONS AND CONCLUSIONS 



The darkon model can have significant effects on collider physics. Since the darkon interacts 
primarily with Higgs bosons, its greatest impact is on the Higgs sector. In either the SM+D or the 
THDM II+D, the existence of the darkon can give rise to huge enhancement of the Higgs width via 
the additional process h —> DD because the darkon-Higgs coupling can be sufficiently large, as can 
be seen from Fig. [3] or |5j In particular, our parameter choice above for the THDM II+D leads to 
the h DD partial width and branching ratio shown in Figs. [7] and El We see that for small m D 
values, the Higgs width is increased considerably and dominated by the h — > DD mode, especially 
if m h < 2m w , in which case the Higgs partial width into SM particles is small. Since the darkon 
is stable, the darkon pairs produced in h — > DD will be invisible, and so in this case the invisible 
branching fraction of h can be substantially bigger than its branching fraction to SM particles. For 
large m D values, the h — > DD contribution becomes small and the h decay is like that in the case 
without the darkon. 

Although the increase of the invisible widths of Higgs bosons will obviously affect Higgs searches 
at the LHC 0,0,8, it will still be able to provide an environment for probing the darkon model. 
For a Higgs boson with a large invisible branching fraction (> 60%) and a mass within the range 
120 GeV < m h ^ 300 GeV, direct searches at the Compact Muon Solenoid (CMS) experiment 
through the usual SM modes may be unfeasible with 30fb _1 of integrated luminosity [8]. However, 
such a Higgs boson can be observed at ATLAS with the same integrated luminosity, via weak-boson 
fusion or Z-Higgstrahlung, by looking for missing energy from the decay jsl, Eo| . Considering one 
of our THDMII+D examples, with tana = 0.015 and tan/3 = 30, we find that h has SM-like 
couplings to the weak bosons and invisible branching ratios higher than 0.6 in much of the m D 
range, as can be seen from Fig. El which serves to illustrate the testability of the darkon model at 
the LHC. Moreover, in the SM+D, if m h > 2m z , say m h = 300 GeV, the total Higgs width can 
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FIG. 7: Partial width of h — > DD in the THDMII+D as a function of darkon mass m D for Higgs mass 
values m h = 120, 200, 350 GeV and the two cases in Fig. [5j Only regions corresponding to X h < 1 have 
been plotted. 
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FIG. 8: Branching ratio of h — > DD in the THDMII+D as a function of darkon mass m D , corresponding 
to T(h -> DD) in Fig. [3 



be measured with a precision of up to 10% at ATLAS with 300 fb" 1 integrated luminosity 2l|, and 
then the darkon contribution can be inferred after a comparison with the SM prediction for the 
width. 

In conclusion, we have studied the experimental constraints from XENON10 and CDMSII on 
the darkon model of dark matter. In the SM+D, a darkon with mass in the ranges of 10 GeV to 
(50, 70, 75) GeV for Higgs boson masses of (120, 200, 350) GeV is ruled out as a WIMP candidate 
by the experimental restrictions. By extending the SM+D to a two-Higgs-doublet model plus 
a darkon, the THDMII+D, the experimental limits can be circumvented due to suppression of the 
darkon-nucleon elastic cross-section at some values of tana tan/3. Future DM search experiments 
can further constrain the parameter space of the model. Using the darkon-Higgs coupling extracted 
from the DM relic density, we have found that the total decay width of h in either the SM+D or the 
THDMII+D can be greatly enhanced by a large contribution from the invisible mode h — > DD. 
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The substantial increase of the invisible decay width of the Higgs boson would lead to a sizable 
reduction of its branching fraction to SM particles. Although this could significantly affect Higgs 
searches at the LHC, we expect that it will still be able to probe the darkon model. 
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search experiment by the TEXONO Collaboration [24| which reported an upper limit on the spin- 
independent WIMP-nucleon elastic cross section for WIMP-mass values below lOGeV, to as low 
as 3 GeV. We have subsequently applied our analysis to darkon-mass values from 3 to 10 GeV and 
found that the two models we considered, the SM+D and the THDM II+D, satisfy the bounds from 
TEXONO. However, their upcoming ultralow energy germanium detector (ULEGe) experiment may 
be sensitive enough to place constraints on the models. 



Appendix A: Derivation of Higgs-nucleon coupling 



In the t-channel Higgs-mediated darkon-nucleon scattering, the energy exchanged is small com- 
pared to the nucleon mass. Therefore, we can use chiral perturbation theory at leading order to 
derive the Higgs-nucleon coupling. 

The relevant Lagrangian is written down in terms of the lightest meson- and baryon-octet fields, 
collected into 3x3 matrices B and £. At leading order in the m s expansion, the chiral strong 
Lagrangian can be written as 



22 



£ s = -m B {BB) + b D (B {M + , B}) + b F (B [M+, B}) + b (M+> {BB) 

+ lf 2 B (M + ) + ••• , (Al) 

where only the relevant terms are shown, m B is the baryon mass in the chiral limit, (• • • } = 
Tr(- • • ) in flavor-SU(3) space, / is the pion decay constant, M + = &M q & + f Mjf with M q = 
diag(m n , m d , m s ) , and b DF0 and B are free parameters which can be fixed from data. In the 
following, we will take the isospin-symmetric limit m u = m d = rh. 

From flAip . one can derive the nucleon mass m N and the pion- nucleon sigma term o~^ N , 
m N = m B - 2(b D + b F + 2b Q )rh - 2(b D - b F + b )m s , 

(A2) 

a nN = rhdm N /dm = —2(b D + b F + 2b )m, 

as well as expressions for the masses of the other hadrons (E, S, K, w) in terms of b DF0 and B . 
Empirically, a nN is not very precisely determined and hence will be a source of uncertainty in our 
calculation. For definiteness, we adopt cr nN = 45 MeV (l7[ . From the expressions for the masses 
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and cr ffJV , one can obtain 



B 



-0", 



2m 2 K + ml (m s + m s ) (2m^ - m^) -2m N m\ 
N 2ml 4(m 2 K -ml) 



(A3) 



where numerically m x , m n , m H , m s , and are isospin- averaged values. 

We can now derive the light-quark contribution to g^NU m ©• Since the terms containing 
M in Eq. (IAip are the chiral realization of the light-quark mass terms in the quark Lagrangian 



- —m u uu — m d dd — m s ss + ■ ■ ■ , we have 

dC dC dC 

(N\k u m u uu + k d m d dd + k d m s ss\N) = -(N\k u m u —^ + k d m d —^ + k d m s —^\N) 

om„ am, am a 



\ (\i k d ) a nN + k d ( 



m N — m 



B 



NN 



(A4) 



using Eq. flA2j) . where the last line has been obtained by averaging over the p and n matrix elements. 



For the heavy-quark contribution to g NNn , we use the relation 14 

m Q QQ = 



_ /~1 (~~<liv 

-3 g7r a^v a 



(A5) 



at leading order in the heavy quark expansion, for Q = c, b, or t, and the nucleon matrix element 
of the trace of the energy-momentum tensor 231 



m N NN = (N\9£\N) = (N\m u uu + m d dd + m s 



9<x 



ss 



G ailu G^\N) 



8tt afiu 



(A6) 



where G% u is the field strength tensor of the gluons. Setting k u = k d = 1 in Eq. (1A4I) and comparing 
it with Eq. (1A6I) . we conclude that 



(iV|m Q QQ|iV) 



2 Qrv 2 - 

■— <iV| — G OUI / GTW = —mo AT AT. 
27 871" 27 



(A7) 



It follows that 



(N\k u m c cc + k d m b bb + k u m t it\N) = — (2k u + k d ) m B N N . 



(A8) 



Combining Eqs. ( 1A4|) and ( 1A8[) . we finally find 



(A9) 



We remark that, for k u and fc d being comparable in size, g pp u and g nnH differ from each other and 
from the isospin- averaged g^Nu above by at most a few percent. 
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